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Abstract 

We provide a summary of new developments in the area of direct reaction the¬ 
ory with a particular focus on one-nucleon transfer reactions. We provide a 
status of the methods available for describing (d,p) reactions. We discuss the 
effects of nonlocality in the optical potential in transfer reactions. The results 
of a purely phenomenological potential and the optical potential obtained from 
the dispersive optical model are compared; both point toward the importance 
of including nonlocality in transfer reactions explicitly. Given the large ambi¬ 
guities associated with optical potentials, we discuss some new developments 
toward the quantification of this uncertainty. We conclude with some general 
comments and a brief account of new advances that are in the pipeline. 

1 Introduction 

Nuclear reactions are an important and versatile tool to study nuclei, particularly those nuclei away from 
stability. One-nucleon transfer reactions can provide information about the single-particle structure of 
the nucleus of interest (see e.g., |[^|^). Deuterons are often used as a probe. However, due to its loosely 
bound nature, it is understood that deuteron breakup effects need to be considered carefully. 

A number of experimental programs in rare-isotope facilities worldwide are also using (d,p) and 
(d,n) reactions as a tool to extract capture rates of astrophysical relevance (e.g., (^Q). Given the many 
ongoing and planned experimental efforts in this direction, it is critical that the reaction theory used to 
interpret those results be reliable and that the uncertainties in the reaction theory be well understood. 

2 Status of the treatment of reaction dynamics 

The treatment of deuteron induced reactions involving intermediate and heavy nuclei pose severe chal¬ 
lenges to microscopic theory. For most reactions involving these systems, the problem is cast as a 
three-body scattering problem with effective nucleon-nucleus interactions. We first briefly review fhe 
currenf sfafus of reacfion fheories for describing deuferon-induced one-nucleon Iransfer reacfions. 

The continuum discrefized coupled channel mefhod (CDCC) Q developed in fhe eighties has 
now been benchmarked againsf fhe exacf Faddeev mefhod Q. The resulfs of fhis benchmark |7| demon- 
sfrafe fhaf while fhe reduction fo one Jacobi componenf is adequafe for Iransfer reacfion af low energies, 
fhe convergence rale for fhe low-energy breakup disfribufions is exfremely slow and does nol allow for 
reliable exfrapolalions. Reacfion calculafions af higher energies inlroduce fhe issue of fhe energy depen¬ 
dence in fhe optical pofenfials. While in coordinafe-based mefhods fhere is fypically a fixed choice of 
the energy at which the optical potentials are evaluated, in the current implementation of the Faddeev 
method Q, calculations are performed in momentum space and can take into account the explicit energy 
dependence of the interaction in the reaction dynamics. This poses an ambiguity in the comparison that 
is discussed in detail in Q . 


While both the CDCC and Faddeev methods are computationally intensive, the adiabatic distorted 
wave approximation (ADWA) f8], which is essentially a simplification of CDCC to make all the channels 
in the continuum degenerate with the ground state (adiabatic approximation), provides a very efficient 
tool to analyze transfer reactions and explicitly includes deuteron breakup to all orders. A benchmark of 
this approach with the exact Faddeev has also been performed and the results show that the adiabatic 
method, while only valid for deuteron-induced transfer reactions, provides as good an agreement as the 
CDCC method. 


Understanding the regions of validity of these reaction theories is very important. However, the 
current implementation of the Faddeev method Q also has its limitations, a fact that became clear during 
the benchmarking process pO} . As the effects of the Coulomb force increase, namely as the energy 
decreases and/or the charge increases, the method of Coulomb screening used by Q begins to fail. 
For example, currently we are not able to produce exact cross sections for ^°Be(d,p) below « 5 MeV, 
^®Ca(d,p) below ss 15 MeV, or any reaction cross section of interest for isotopes heavier than Ni. One 
way to overcome this problem is to cast the Faddeev equations in the Coulomb basis instead of the 
traditional plane wave basis O- By doing this, one avoids Coulomb screening completely. The TORUS 
collaboration a topical collaboration in nuclear theory with the goal of advancing methods of (d,p) 
reactions, has put in place many pieces of the puzzle necessary to solve the Faddeev equations in the 
Coulomb momentum space basis 113 -151. We expect in the next couple of years to have a completely 
new and fully developed code that can address the three-body problem of A(d,p)B for nuclei with large 
Z. 


3 Nonlocality in optical potentials and transfer reactions 

One of the most important ingredients in predictions of reaction cross sections are the optical potentials. 
From a microscopic point of view, it is clear that the nucleon-nucleus optical potential should be non local 
due to many-body effects, including antisymmetrization and excitations. However, phenomenological 
potentials have traditionally been assumed local, for convenience. The consequence of this simplification 
is that it becomes strongly energy dependent. As many-body methods improve, including more and more 
correlations and a careful treatment of the continuum, one can expect the optical potential to be fully 
extracted from these many-body theories in the near future. It is therefore timely for reaction theory to 
be prepared to handle nonlocality in the effective interactions and to understand the magnitude of the 
effects. 


It is clear that for any non-local potential, one can construct local phase-equivalent potentials, i.e., 
reproducing the exact same elastic scattering. Even when the elastic scattering is exactly reproduced, 
nonlocality can imprint itself in other channels. In | |T^ 171 the problem is addressed for one-nucleon 
transfer reactions. In the first of these works 1161, a large number of reactions at different beam energies 
and for different targets are considered. The Percy and Buck non-local optical potential was used 1181 and 
local phase-equivalent potentials were constructed. Finally, the inclusion of nonlocality in one-nucleon 
transfer reactions populating single-particle states was studied. The study was performed in DWBA and 
nonlocality was included only in the exit channel, namely in the neutron-bound state and the proton- 
distorted wave. A similar study was performed using the dispersive optical model | [T^ . Generally, the 
nonlocality in the scattering state reduced the cross section, while the nonlocality in the bound state 
increased it, with the net effect being an increase as large as 30%. In we consider hole states in 
40Ca, again using DWBA and nonlocality only in the exit channel. The dispersive optical model and the 
Percy and Buck potentials provide similar results. The effect of nonlocality on transfer is very large for 
these hole states (up to 50%). 


In Fig. I we show the predicted transfer cross sections for ^^®Sn(d,p)^^^Sn at = 20 MeV. 
Comparing the local DWBA predictions (blue dotted line) with the corresponding non local (green dot- 
dashed line), we see an effect of up to 20% in magnitude. Also shown, for completeness, are the same 
results within ADWA (red dashed for the local and black solid for the non local). While ADWA and 
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Fig. 1: Angular distributions for ^^®Sn(d,p)^^’^Sn at Ed = 20 MeV; comparing the inclusions of nonlocality in the 
exit channel for ADWA and for DWBA. 



Fig. 2: Angular distributions for ^^®Sn(d,n)^^^Sb at = 20 MeV; comparing the inclusions of nonlocality in the 
exit channel for ADWA and for DWBA. 
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Fig. 3: Contour plot for the residuals resulting for the fit of neutron elastic scattering on "^®Ca at 23 MeV, when 
varying the radius parameter r of the real part of the optical potential and the corresponding depth V. 


DWBA show noticeably different angular distributions (a difference coming from the different tteatment 
of deuteron breakup), the relative effect of nonlocality is magnified in ADWA. Note that in these ADWA 
calculations, the nonlocality in the deuteron channel is not yet included. Nonlocality in the deuteron 
channel will be discussed elsewhere | [20| ]. 

Because nonlocality is introduced in the nuclear interaction only, one would expect to find similar 
effects in the study of (d,n) reactions. In Fig. 2 we show the angular distributions for ^^®Sn(d,n)^^^Sb at 
Ed = 20 MeV. All lines in Fig. 2 correspond to DWBA calculations. The result when using the local 
interactions is shown by the red dashed line and that including the non-local interactions is the black 
solid line. Note that the ground state of is different than that of ^^^Sn, given the different proton 
and neutron numbers of the target. While the ground state of ^^^Sn is a lfiii /2 state, in the ground 
state is lgj /2 with a significant difference in the binding energy. In addition, the proton state is subject 
to Coulomb repulsion. All these differences compound to produce a very different effect of nonlocality 
for the (d,p) and (d,n) reactions on heavy nuclei. For ^^®Sn(d,n)^^^Sb, the effect of nonlocality is very 
large, doubling the cross section at the peak. We also show in Fig. 2 the separate effects of introducing 
nonlocality in the neutron bound state (green dot-dashed line) and the proton scattering state (blue dotted 
line). Nonlocality in the proton scattering state reduces the cross section but only by a little. The larger 
effect is in the nonlocality in the bound state. 

These studies all call for the need to include nonlocality explicitly in reactions for a reliable de¬ 
scription of the process. 



Depth [MeV] 

Radius [fm] 

Diffuseness [fm] 

Volume Term (real) 

162.4 (103.2) 

0.98 (1.05) 

0.69 (0.86) 

Surface Term (imaginary) 

48.7 (16.5) 

1.1 (1.43) 

0.30 (0.67) 

Spin-Orbit Term 

3.33 

1.07 

0.66 


Table 1: Best fit parameters for the elastic scattering of ^®Ca(d,d)'*®Ca. The spin-orbit term was not fit. All optical 
model potential parameters not listed here were zero. The optimization was initialized at the value in parentheses. 
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Fig. 4: Elastic scattering angular distribution for neutrons on "^®Ca at 23.2 MeV: solid black line is the best fit, the 
brown corresponds to the 95% confidence band and in black are the data points 


4 Optical potential uncertainties 

The study in Q demonstrates that the optical potential uncertainties are reduced in ADWA when com¬ 
pared to DWBA, the main reason being that ADWA relies only on nucleon optical potentials, whereas 
DWBA requires a deuteron optical potential that is less well constrained. Knowing that these effective 
interactions carry ambiguities, it is important to quantify the uncertainties in the predicted reaction ob¬ 
servables. So far estimates of these uncertainties have been performed by comparing the results when 
two arbitrary optical potentials are chosen (e.g., Q). Then the results of the estimated error bars depend 
on the choice of the two representative potentials. Alternatively, one can make use of modern techniques 
in the field of uncertainty quantification and explore these tools in the domain of reaction theory. 


Our plan is to use modern statistical tools to quantify uncertainties in nuclear reactions (see 
e.g., I^). This is largely an untapped field and therefore we start with a simple case. We look at 
the case of ‘^^Ca(d,p)^®Ca in DWBA and examine the effect of the uncertainty from constraining the 
deuteron optical potential on the predicted transfer cross section. We use elastic scattering data at 23.2 
MeV from |23| and perform a minimization to find the best fit. We take a typical Woods-Saxon 
shape and the parameterization of 1241 as a starting point for the minimization routine. We allow all 
three parameters in the central volume real and the central surface imaginary terms (depth, radius, and 
diffuseness) to vary simultaneously (the spin-orbit interaction is kept fixed). Parameters are only allowed 
to vary within physical limits. The minimum found is summarized in Table I, with initial parameters 
shown in parentheses. We then verify that the distribution of the total in the elastic channel, around 
the minimum, is Gaussian. This can be seen by the elliptical behavior of the | |x^| I contour plots. In Fig. 

we show an example of such a plot, where the vertical axis is the potential depth and the horizontal 
axis is the radius of the volume term. 


We then pull 200 sets randomly from the Gaussian distribution based on the elliptical behaviour 
surrounding the minimum, throw out the largest 5 and lowest 5 predicted cross sections per angle, to 
obtain the 95% confidence band. Using the same 200 sets, we define the 95% confidence bands indepen¬ 
dently for elastic and transfer. The results are depicted in Figs. |^an d0 

The error bars considered in quantifying the x^ for these calculations were the experimental error 
bars. However, often the difference between data and theory is larger than the experimental error bars. 
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Fig. 5: Angular distribution for ^®Ca(d,p)‘*®Ca at 23.2 MeV: solid black line is the best fit, the brown corresponds 
to the 95% confidence band. 


One can then repeat the procedure including the theoretical error. Also noteworthy, there are cases in 
which the distribution around the minimum is not Gaussian. Then one may need to pull the random sets 
from the regions of constant directly and/or appeal to Bayesian techniques |25 261. 


A refinement of the fitting can be done by including a larger set of data, with perhaps a variety 
of observables. As an example, we have started to look at simultaneously fitting elastic and inelastic 
cross sections, including target excitation in the reaction model. While a larger number of parameters are 
available for fitting, the actual procedure for extracting confidence bands is fhe same. A more challenging 
problem is fhaf of fhe uncerlainfy in fhe fheorefical model ilself. More work will be needed fo defermine 
fhe besf approach. 


5 Outlook 

As described here, fhere are a number of excifing developmenfs in reaction fheory. The implemenfafion 
of fhe Faddeev equafions in fhe Coulomb basis is progressing and will provide reliable predictions for 
(d,p) reacfions for a wide range of nuclei and energies. The sfudies performed on non-local inferacfions 
demonsfrafe fhe need fo include fhese inferacfions explicifly in describing fransfer processes. While fhese 
are crifical sfudies for fhe fufure of fhe field, fhey still rely on fhe knowledge of fhe effective inferacfions 
befween nucleons and fhe large!. 

Allhough Iradifionally fhese effeclive inferacfions have been exfracfed phenomenologically, we 
are currenfly beginning an era where fhe level of nuclear many-body mefhods may enable fhe exlraclion 
of fhese pofenlials direclly from fhe NN inleraclion. A new collaboration befween our group and fhe 
ORNL fheory group aims al exlracling fhese pofenlials from coupled-clusler fheory. While one mighl 
slill expecl a deficiency in fhe absorptive componenl, fhe inclusion of a variety of correlafions and fhe 
confinuum basis hold promise for using Ibis melhod versus olher many-body mefhods. 

Lasl bul nol leasl, we are sfarling fo develop fechniques fo quanlify uncertainties in reaction fheory. 
This work is still in ils infancy. As described in Ihis proceedings, we are still exploring fhe besf way fo 
characlerize fhe uncerfainfies for fhe limifed problem of fhe effecl of fhe errors in a given dafa sel in 
conslraining fhe inleraclion and fhe propagation fo reaclion observables. There are many ofher sources 
of uncerfainfies in reaclion fheory, which may be equally or even more imporlanl. These need fo be 
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considered case by case and then accounted for jointly so that we can provide a reliable and useful 
procedure for uncertainty quantification in this field. 

This work was supporfed by fhe Nafional Science Foundafion under Granf No. PHY-1403906 
and PHY-1304242 and fhe Deparfmenf of Energy, Office of Science, Office of Nuclear Physics under 
award No. DE-FG52-08NA28552, de-sc0004087, DE-87ER-40316. This work was also supporfed by 
fhe U.S. Deparfmenf of Energy, Office of Science, Advanced Scienfific Compufing Research program 
under confracf number DE-AC02-06CH11357. 
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